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Conversion  Factors,  Non-SI  to  SI  (Metric) 
Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

feet 

inches 

kilotons 

pounds  (force)  per 
square  inch 

pounds  (mass) 

pounds  (mass)  per  cubic 
foot 

pounds  (mass)  per  foot 

square  feet 

tons  (mass)  per  cubic 
foot 

tons  (mass)  per  square 
foot 


By 

0.3048 

25.4 

4.184 

0.006894757 

0.4535924 

16.01846 


_ To  Obtain _ 

metres 

millimetres 

terajoules 

megapascals 

kilograms 

kilograms  per  cubic  metre 


1.488164 

0.09290304 

32,036.978 


kilograms  per  metre 
square  metres 
kilograms  per  cubic  metre 


9,764.856 


kilograms  per  square  metre 
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IMPORTANCE  OF  . RIGID-BODY  MOTIONS  ON  A 


SINGLE-DEGREE-OF-FREEDOM  MODEL 


Purpose  and  Scope 

1.  The  purpose  of  this  study  is  to  determine  the  effect  of  rigid-body 
motion  on  the  single-degree-of-freedom  (SDOF)  model  presently  used  for  dynamic 
analysis  of  buried  roof  slabs.  The  effect  of  rigid-body  motion  is  incorpo¬ 
rated  by  adding  a  second  degree  of  freedom  to  the  model.  The  second  degree 
will  model  the  motion  of  the  entire  structure,  which  is  controlled  by  the 
dynamic  pressure  applied  and  the  resistance  of  the  soil  mass.  The  results  of 
this  model  should  not  be  interpreted  to  be  more  accurate  than  those  produced 
by  the  SDOF  model.  This  model  is  valid  only  for  studying  the  effect  of  rigid- 
body  motion  since  this  model  has  no  experimental  data  or  experience  to  support 
its  results. 

2.  The  results  that  will  be  examined  are  the  relative  displacements 
between  the  roof  and  the  structure.  The  absolute  roof  displacements  will 
become  larger  as  a  result  of  the  rigid-body  motion.  However,  the  relative 
displacements  between  the  roof  and  structure  are  the  only  displacements  that 
result  in  damage  to  the  roof.  Variations  in  soil  resistance  and  loading 
intensity  are  examined. 


Procedure  Using  Two-Degrees-of-Freedom  Model 


3.  The  two-degrees-of-freedom  model  to  be  used  for  this  study  are  shown 
in  Figure  1 . 

l  H(U mi 

R2  (SOIL  RESISTANCE) 

fA^j  (SOIL  MASS)  ____ 

F2  <t)  f  V2 

Rt  (ROOF  RESISTANCE) 


[M  j  (ROOF  MASS) 

Fi(t) 


~F 

Figure  1.  Degree-of-freedom  models 
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4.  The  roof  resistance  is  a  bilinear  function  based  on  previous  SDOF 
models.  This  bilinear  function  is  shown  in  Figure  2. 


y2  -  vi  . in. 

Figure  2.  Roof  resistance 


5.  The  loading  functions  F^  and  F2  are  dependent  on  the  pressure¬ 
time  history  of  the  weapon.  The  pressure-time  (P(t))  relation  for  the  first 
50  msec  is  given  in  Figure  3.  The  total  pressure  duration  is  taken  to  be 
250  msec,  and  the  pressure  is  assumed  to  decay  linearly  to  zero  from  50  to 
250  msec.  Based  on  the  area  of  roof  and  foundation,  the  forcing  functions  are 
calculated  as  follows: 

F 1  ( t )  =  P(t)  (109.5M240M0.66)  =  17,344.8  P(t)  1b* 

(dimension  of  roof)  (soil  arching) 

F2(t)  =  P(t)( 1 45 . 5 ) ( 264 )  =  38,412  P(t)  lb 
(assumed  foundation  dimensions) 

Soil  arching  is  the  ability  of  a  soil  to  transfer  loads  from  one  location  to 
another  in  response  to  a  relative  displacement  between  the  locations.**  The 

p 

total  mass  of  the  roof  is  855  lb-sec  /ft.  The  minimum  amount  of  contributing 

soil  mass  from  the  foundation  was  assumed  to  be  5  ft  below  the  foundation. 

This  results  in  the  mass  of  the  foundation  being  approximately  7,000  lb- 
2 

sec  /ft.  This  mass  is  based  on  assuming  that  the  foundation  material  has  a 
weight  of  110  lb/ft^.  The  foundation  resistance  is  based  on  Terzaghi's  co¬ 
efficient  of  vertical  subgrade  reaction. ^  Since  a  cohesive  material  is  used 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 

**  3.  A.  Kiger,  T.  R.  Slawson,  and  D.  W.  Hyde.  1984  (Sep).  "Vulnerability 

of  Shallow-Buried  Flat-Roof  Structures;  Final  Report:  A  Computational  Pro¬ 
cedure,"  Technical  Report  SL-80-7,  Report  6,  US  Army  Engineer  Waterways  Ex¬ 
periment  Station,  Vicksburg,  Miss. 

f  Karl  Terzaghi.  1985.  "Evaluation  of  Coefficients  of  Subgrade  Reactions," 
Geotechnique,  Vol  5,  pp  297-326. 
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around  the  structure,  values  for  stiff  to  very  stiff  clay  were  used  to  develop 
the  bilinear  resistance  curves  for  the  foundation,  the  calculations  for  these 
curves  are  shown  in  Appendix  A. 

6.  Two  different  limiting  bilinear  curves  can  be  produced  using  the 
summary  of  the  data  from  Appendix  A  (Figure  4).  The  and  are 

based  on  assuming  the  foundation  material  to  behave  between  a  stiff  to  very 
stiff  clay.  These  are  realistic  values  that  will  limit  the  valid  ranges  of 
the  soil  stiffness.  However,  according  to  Terzaghi,  these  values  are  used  for 
loads  not  to  exceed  one-half  of  ultimate  bearing  capacity  pressure  of  the 

F  =  13  x  106  LB 

(FOR  MAX  SURFACE 
PRESSURE  =  600  PSD 

Kmax  =  6.7x106  LB/IN. 

Kmin  =  3.4  x  106  LB/IN. 


v2,in. 


Figure  4.  Soil  resistance 
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structure  and  also  for  static  loads  applied  only  to  structure.*  The  ultimate 
bearing  capacity  is  increased  with  the  addition  of  a  surcharge  load  and  thus 
makes  the  soil  resistance  dependent  on  time.  Because  the  response  of  the 
foundation  is  due  to  dynamic  forces  and  the  force  is  of  short  duration,  the 
maximum  pressure  was  assumed  for  the  purpose  of  calculating  the  ultimate 
bearing  capacity  to  always  be  present  at  the  surface. 

7.  In  addition  to  changing  the  properties  of  the  soil  resistance,  the 
effect  of  the  rigid-body  motion  will  also  be  evaluated  for  additional  soil 
mass  and  increased  dynamic  loading.  Therefore,  load  cases  3  and  4  in  Table  1 
are  identical  to  load  cases  1  and  2,  except  for  the  increased  soil  mass. 

According  to  Clough  and  Penzien,**  the  added  mass  for  a  circular  section  is 

3 

1.5  pr  ,  where  p  =  weight  per  unit  volume  and  r  =  radius  .  The  mass 

3 

for  load  cases  3  and  4  is  1.5  pb  ,  where  b  =  the  minimum  foundation 
dimension  . 

8.  The  relative  magnitude  of  the  forcing  function  is  also  a  major  vari¬ 
able  controlling  the  dynamic  response  of  a  structure.  Therefore,  load  cases  5 
through  8  are  identical  to  load  cases  1  through  4,  except  that  the  load  func¬ 
tions  F^t)  and  F2(t)  have  been  multiplied  by  a  factor  of  3  and  the 
maximum  force  for  the  soil  resistance  curve  is  36  x  10^  lb  (see  Appendix  A). 

9.  Load  cases  9  and  10  have  large  masses  for  the  foundation  and  stiff 
soil  resistance  curves.  These  factors  will  effectively  eliminate  any  rigid- 
body  motion  and  will  produce  displacements  to  compare  with  other  load  cases. 
Load  cases  9, should  be  compared  with  cases  1  through  4  to  determine  the  effect 
of  rigid-body  motion  under  relatively  low  pressures  and  light  damage  predic¬ 
tions.  The  results  from  load  case  10  can  be  used  to  evaluate  load  cases  5 
through  8,  which  had  relatively  high  overpressures  and  heavy  damage  predic¬ 
tions.  Table  1  gives  all  of  the  variables  for  each  load  case  considered. 

10.  Table  2  summarizes  the  results  of  dynamic  analysis  of  the  two- 
degrees-of-freedom  models.  Plots  displaying  the  displacements  of  the  struc¬ 
ture,  foundation,  and  the  relative  displacements  are  shown  in  Appendix  B. 

11.  As  seen  from  Table  2,  the  rigid-body  motion  of  the  foundation  and 
structure  can  significantly  affect  the  relative  displacements  between  the 
structure  and  the  roof.  When  load  cases  1  and  2  are  compared  with  load 

*  Karl  Terzaghi  and  Ralph  Peck.  1967.  Soil  Mechanics  in  Engineering  Prac¬ 
tice  ,  Wiley,  Mew  York. 

**  Ray  W.  Clough  and  Joseph  Penzien.  1975.  Dynamics  of  Structures,  McGraw- 
Hill,  Mew  York. 
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Table  1 

Model  Parameters 


Load 

Case 

Mi 

lb-sec2/ 

in. 

m2 

lb-sec2/ 

in. 

RIM  Max 
Resist¬ 
ance 
106  lb 

fiK1 

10°  lb/ 
in. 

R2M  Max 
Resistance 
Soil,  106  lb 

k2 

104  lb/ 
in . 

Maximum 

Pressure 

psi 

1 

71.25 

588 

5.8 

12 

13 

2.9 

610 

2 

71.25 

588 

5.8 

12 

13 

5.8 

610 

3 

71.25 

1,200 

5.8 

12 

13 

2.9 

610 

4 

71.25 

1,200 

5.8 

12 

13 

5.8 

610 

5 

71.25 

588 

5.8 

12 

36 

2.9 

1,830 

6 

71.25 

588 

5.8 

12 

36 

5.8 

1,830 

7 

71.25 

1,200 

5.8 

12 

36 

2.9 

1,830 

8 

71.25 

1,200 

5.8 

12 

36 

5.8 

1,830 

9 

71.25 

3,000 

5.8 

12 

9,000 

9,000 

610 

10 

71.25 

3,000 

5.8 

12 

9,000 

9,000 

1,830 

Table  2 

Results  of  Dynamic  Analysis  of  the  Two-Degrees-of-Freedom  Models 


Load 

V  "2 

in. 

Max  Pressure 

Case 

Soil  Type 

psi 

1 

1.04 

Stiff 

610 

2 

1.01 

Very  stiff 

610 

3 

1.57 

Stiff 

610 

4 

1.59 

Very  stiff 

610 

5 

118.0 

Stiff 

1,830 

6 

114.0 

Very  stiff 

1,830 

7 

108.0 

Stiff 

1,830 

8 

113.0 

Very  stiff 

1,830 

9 

2.66 

Stiff 

610 

10 

109.0 

Very  stiff 

1,830 

case  9, 

the  maximum  percent  difference 

is  89  percent. 

When  the  mass  is  in- 

creased  to  a  more  reasonable  value  (as  in  load  cases  3  and  4),  the  percent 


difference  with  load  case  9  becomes  51  percent.  However,  when  the  loads  are 
increased,  the  percent  difference  decreases,  as  seen  when  load  cases  5  and  6 
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are  compared  with  load  case  10.  This  comparison  results  in  a  7.9  percent  dif¬ 
ference  between  load  cases  7  and  10.  When  the  mass  is  increased,  the  percent 
difference  becomes  1  percent  between  load  cases  7  and  10.  The  two  principal 
factors  controlling  the  relative  displacement  are  the  forcing  function  and  the 
assumed  amount  of  foundation  that  acts  with  the  structure. 

Conclusions 

12.  These  results  indicate  that  at  overpressures  required  to  cause  se¬ 
vere  damage  to  a  relatively  hard  shallow-buried  structure,  rigid-body  motion 
will  have  little  effect  on  the  damage  predicted  from  an  SDOF  analysis.  How¬ 
ever,  at  relatively  low  overpressures,  rigid-body  motion  can  result  in  sig¬ 
nificantly  less  damage  than  would  be  predicted  by  an  SDOF  analysis.  Since 
typical  design  calculations  are  for  relatively  low  overpressures  (i.e.,  light 
damage),  it  is  important  to  include  the  effects  of  rigid-body  motion.  How¬ 
ever,  for  vulnerability,  calculations  that  predict  pressure  required  to  fail 
the  structure,  including  the  effects  of  rigid-body  motion,  may  not  be 
necessary. 


9 


Calculations  for  Bilinear  Soil  Resistance  Curves 


Values  for  Subgrade  Modulus* 


Terzaghi's 
Consistency  of  clay 
Ec 


Stiff  Very  stiff 

75  tons/ft^  150  tons/ft^ 


Values  valid  for  pressures  <  one-half  of  the  ultimate  bearing  capacity  qd 
where : 

qd  =  1.2  cNc  +  (yDf  +  q)  Nq  +  0.4  y  BN^** 
y  =  Unit  weight  of  soil  =  115  1b/ft^ 

Df  =  Height  of  overburden  =  10  ft 
B  =  Width  of  foundation  =  12  ft 
C  =  Cohesion  =  1/2  qu 
q  =  Surcharge  load  =  600  psi 


qu  =  Unconfined  compressive  strength 
Consistency  of  clay  Stiff 

qu  (tons/ft^)  1-2 

C  (tons/ft^)  0.75 

For  0=0 

Nc  =  {2  +  it) 

V-  1 

N  =  0 


Very  stiff 
2-4** 
1.5 


Qd, 


min 


qdmin 

=  1.2  (0.75)  (5.14) 

(2,000)/144 

qdmin 

=  64.27  +  607.99 

q^min 

=  672.26  psi 

qdmax 

=  1.2  (1.5)  (5.14) 

(2,000)/144 

=  128.5  +  607.99 

qdmax 

=  736.49  psi 

*  qdmax 

=  700  psi 

n^.QQ)  (600) 


144 


(115)  (10) 


144 


(600) 


Note:  Bearing  capacity  is  controlled  by  the  loading  on  ground  surface. 


Maximum  allowable  force  _  qDA  _  700  ( 1 45 . 5 ) ( 264 )  _  n  ik 

[p<t>max  --  6°°  Psi]  '  2  ‘  2  * 


*  Karl  Terzaghi  and  Ralph  Peck.  1967.  Soil  Mechanics  in  Engineering  Prac¬ 
tice,  Wiley,  New  York. 

**  Ray  W.  Clough  and  Joseph  Penzien.  1975.  Dynamics  of  Structures,  McGraw- 
Hill,  New  York.  . ‘ 


A2 


Maximum  allowable  force  _  ( 1 ,900)(  145 .5) (264)  _  r,a  m6  lu 
[P(t)maX  =  1,800  psi]  -  2  -  3b  x  10  lb 

Spring  stiffness,  K  lb/in. 


K  = 


K  = 


E  (connection  for  rectangular  footings)  area 

3 _  _ _ 

(connection  for  larger  than  1x1) 

E  (m  +  0.5)  A 
s 


B  1.5  m 

M  =  ratio  of  rectangular  footing  sides 


Kmax  = 


’(max) 


(m  +  0.5)  A 


B  1.5  m 


.  ( 150) (2, 000) (0.551 1  +  0.5)  m5  5U26i|) 
"  (1,728)(1, 455)1. 5(0. 5511) 


(min) 


(m  +  0.5)  A 


in 


B  1 .5  m 


=  5.8  x  104  lb/in. 

(75)  (2,000)  (0 .551 1  +  0.5)(  145.5)(26*1) 
( 1 , 728 )  ( 1 145 . 5 )  ( 1  -  5 )  ( 0 . 55 1 1 ) 

=  2.9  x  1041b/in. 
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Appendix  B 

Results  of  Two-Degrees-of-Freedom  Models 


1 


B2 


Figure  B1 .  Displacements  form  two-degrees-of-freedom  model  for  case 


Figure  B2.  Displacements  form  two-degrees-of-freedom  model  for  case 
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Figure  B3.  Displacements  form  two-degrees-of -freedom  model  for  case 


tine,  seconds 


129 


Displacements  form  two-degrees-of-freedom  model  for  case 


Figure  B6.  Displacements  form  two-degrees-of-freedom  model  for 


*>! 
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Figure  B7.  Displacements  form  two-degrees-of-freedom  model  for  case 
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TIME .  SECONDS 


Figure  B9.  Displacements  form  two-degrees-of-freedom  model  for  case 


i2* 


Figure  BIO.  Displacements  form  two-degrees-of-freedom  model  for  case  10 


